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Abstract 

A discrete in time model of ideological competition is formulated tak- 
ing into account population migration. The model is based on interactions 
between global populations of non-believers and followers of different ide- 
ologies. The complex dynamics of the attracting manifolds is investigated. 
Conversion from one ideology to another by means of (i) mass media in- 
fluence and (ii) interpersonal relations is considered. Moreover a different 
birth rate is assumed for different ideologies, the rate being assumed to 
be positive for the reference population, made of initially non-believers. 
Ideological competition can happen in one or several regions in space. 
In the latter case, migration of non-believers and adepts is allowed; this 
leads to an enrichment of the ideological dynamics. Finally, the current 
ideological situation in the Arab countries and China is commented upon 
from the point of view of the presently developed mathematical model. 
The massive forced conversion by Ottoman Turks in the Balkans is briefly 
discussed. 
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1 Nonlinear systems, population dynamics, and 
opinion formation 

Competition among ideologies for the people minds is an interesting phenomenon 
that can lead to large consequences for the social structure and for the social 
evolution of populations. Religious wars in the past or social revolutions during 
the XX-th century are closely connected to such ideological competitions. Ide- 
ological competition is possible because people might and can accept or drop 
ideas and systems of ideas and thereafter act according to the newly adopted 
system of ideas [TJ [2] ■ Population of believers or non-believers do form. The 
number of individuals in these populations can change with time. Thus the 
ideological competition is closely connected to population dynamics state but 
the human ideological struggle has its specific features: animals cannot change 
their kind (zebras cannot become lions; electrons cannot become bosons) but 
people can (though sometimes may not) change their ideology. 

The systems studied as population dynamics are among the best examples of 
nonlinear systems in quantitative social and natural sciences [3| . Different kinds 
of nonlinearities arise in the model equations of such systems, for example (i) 
because of interactions among the individuals or the populations or (ii) because 
of limitations in the environment [4j El [7J [H [9] . Such nonlinear model systems 
require applying methods of nonlinear dynamics [Till [HI EH HI] , chaos theory 
[T51 Hoi [T71 [T5] , theory of stochastic processes [JjJl HOI 121] and even variational 
and other methods of the theory of turbulence 22, 251 l2~i]. 

In this paper, ideological struggle will be considered to occur on a two- 
dimensional surface (regions in a country or regions containing an arbitrary 
number of countries). We shall assume that the population of the correspond- 
ing territory is divided into sub-populations of followers of several ideologies and 
one population of non-believers. Any newborn, whatever each parent ideology is 
necessarily at first a non-believer; any ideology is acquired later. This assump- 
tion has as a consequence that the rate of change of the number of individuals 
by birth and death can be positive only for the population of non-believers. 
The corresponding rate for the sub-populations of the followers of the different 
ideologies is necessarily negative as, for these, only death is possible. 

What will be apparently new in our model is the possibility of migration 
of groups of non-believers or followers of ideologies. Indeed it is going to be 
assumed that non-believers or followers of different ideologies can move from one 
region to another one, even if the regions do not have a common boundary: thus 
allowing motion by airplanes or due to a transit crossing of a region by means 
of train, ship, bus, car, bicycle, horse, camel, donkey, etc.. This movement can 
change its intensity in time; moreover the birth rates and conversion coefficients 
are allowed to change as well. The result of this is a very complicated dynamics 
in time and space. 

The model system of nonlinear equations is taken as a system of difference 
equations which means that one assumes discrete time. One fundamental reason 
for this is our wish that the model be close to the true situation as much as 
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possible, in the sense that data is usually collected for discrete time intervals 
and for a finite number of regions. These two features of the real social life are 
thus accounted for in the model below. 

The paper is organised as follows. In Sect. [2j the system of model equations 
is formulated for the case of a single region, for simplicity and for outlining 
the notations. The possibility of time change of the coefficients of the model 
equation leads us to the notion of special points. Special points of the system of 
these model equations are calculated on the basis of the rules used for calculating 
so called fixed points. However as the coefficients of the system vary in time, 
the points are not truly fixed but move in the phase space. Because of this, 
we use for such points the name special points. If the coefficients of the model 
system of equations do not change in time, the special points become usual fixed 
points. 

Sect. [3] is devoted to a discussion of the model for the case of two, obviously 
adjoining, regions. The possibility of migration being allowed, the dynamics 
become more complicated. Two particular cases are discussed: massive leaving 
of non-believers from one of the regions to the other, in Sect. 3.1 and massive 
invasion by followers of some ideology from one of the two regions to the other 
in Sect. |3.2| Sect. [4] is devoted to the formulation of the model for an arbitrary 
number of regions and an arbitrary number of ideologies. On the basis of the 
general model we discuss the possibilities for the strengthening or weakening of 
an official ideology. Sect. [5] contains several pertinent remarks. Examples are 
taken from recent and non-recent history : Arab countries and China, and the 
spreading of Islam in Europe from the beginning of the I5th century though the 
reader can be aware of other cases according to its own knowledge of history. 
Other discussions and conclusions are found based on a summary in Sect. [6] 



2 Single region 
2.1 Model 

In order not to overload the basic equations with many terms, let us first present 
the model and some discussion in the case of a single region. In this region, con- 
sider that there are n+1 population types, or ideologies: one specific population 
of non-believers, or having no opinion at all about some subject, with spatial 
density p of members at time t, and n populations with spatial densities p\ 
of followers of different ideologies (i — 1, ...,n) about some given question (the 
measure of all densities is people per square kilometer). There is no restriction 
to considering only political ideologies; a religion can be taken as an ideology 
as wellj^] Notice that i can be considered as an index indicating the component 
of the degree of freedom vector characterizing an individual. However, in view of 
being more general, one can consider that some i corresponds to a population, 

Languages are not here considered to be ideologies, due to our consideration that one can 
only belong to one ideology at a given time, thus we exclude multilingual populations for our 
present considerations. 
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having some mere (so called scientific) degree of freedom for characterising part 
of the population, i.e. e.g. characterised by by a political and a philosophical 
attitude, like distinguishing between (i — 1) socialist Buddhists and (i = 2) 
liberal Buddhists, or (i = 3) socialists jews and (i — 4) liberal jews, etcP] Thus 
the total density at time t is 

n 
i=l 

The model (difference) equations for describing the interactions among the 
above n + 1 populations are assumed to be of the Lotka-Volterra type 

n n 

P? +1) - pf = + £ flfpf + £ &gVfVf (2) 

i=o j=o 

where r-*^ is the "biological" rate for changing the density of the corresponding 

population by births and deaths. In other words r\ accounts for the finite 
duration of the life of the people. In the present work, no time delay, memory 
or anticipation is considered between various causes and effects. The change in 
population number occurs at each unit of time t which is discrete. As already 
mentioned in the previous section below we assume that only the population of 
non-believers can have a positive rg rate: thus we assume that (i) all newborn 
babies are non-believers and (ii) for this population, the number of births can 
be larger than the number of deaths. For the other populations of believers, 
the rates rf \ (n > i > 1) are taken to be negative, since only death is possible 
there. 

In addition to births and deaths, we assume in the second and third term of 
the equation r.h.s. that two kinds of conversion exist in the system of popula- 
tions. The first kind of conversion is triggered by the mass media: newspapers, 
Internet, TV, radio, etc. In the model equation for the i-th population, the 
mass media conversion is accounted for by the term fjfpj where /y is the 
corresponding coefficient of conversion influencing the number of followers of 
the j-th ideology who are converted to the i-th ideology because of mass media 
influence. 

The second kind of conversion is the contact conversion that happens, e.g., 
after conversations between members of different ideologies. We assume for 
simplicity that the result from this kind of conversion is proportional to the 
densities of the conflictual contacting ideologies. Thus this kind of conversion 
is modelled by the term fry pf 1 p)p , where bfj is a coefficient that accounts for 
the intensity of the conversion. 

The dynamics of the system of n + 1 populations as described by Eq. ^ can 
be discussed in an n + 1-dimensional Cartesian phase space. The coordinates 
on the axes spanning this space correspond to the densities po, pi, ■ ■ ■ , p n of the 

2 Of course other notations can be used, like Pi t j t k I,... where i indicates the ideology, j the 
religion, k the sex, I the social status, . . . , etc. 
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adepts of the corresponding ideology. The coordinates po, p\ , ■ . . , p n define a 
point in this phase space denoting the social state of the investigated system at 
a given time. As the time t evolves, the system state follows a trajectory: the 
so called system trajectory. This system trajectory can be attracted to different 
kinds of manifolds in the phase space. These manifolds are called attractors: 
they can be points, or more complicated structures. Most of the time the 
system is on its way to some attracting manifold, which in most cases can be 
some equilibrium state. 

2.2 Spreading of ideologies 

Let us now discuss the simplest case, i.e. spreading of one ideology in one 
area due to opinion competition. For simplifying the presentation, let only two 
populations/opinions/ideologies existing: (i) the non-believers po and (ii) the 
followers p\ of ideology 1. The model equations read 

PT X) ~ pf = rfpf + ± fjgpf + ± b^p^pf (3) 

j=0 3=0 



p[ t+1) - pf - r?p? + £ /g>pf + £ b^pfpf (4) 

3=0 j=0 

Without any loss of generality, let also / o = /n = and b 00 = &n = 0, i.e. 
the mass media influence and binary conversion do not increase the number 
of members of both populations per se through some miraculous multiplica- 
tion/addition. Then the system of equations reduces to 

^-tf^tf+ffltf+Wtf (5) 

^-pi'^riV+ZS^+^W- (6) 

Remember that is imposed to be negative and t*q positive. 

The special phase space points connected to the dynamics of the system, Eqs. 
([5| and Q , are obtained when the left-hand sides of the two equations ^ , (|6| 
are equal to 0. We shall call such points special points instead of fixed points 
because of the following reason. Fixed points do not move in the phase space, 
but the special points can move in this space, because the system parameters 
can change with time. Only the special point corresponding to the extinction 
of all populations is a fixed point for all values of the system parameters. 

The (two) special points for the system of equations ([5| and Q are 

(*) (*) n rv\ 
Po = Pi = ( 7 ) 
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and the non-trivial special point 



At) At) _ At) (t) At) At) _ At) <t) 

(t) _ J 10 J 01 ' 1 '0 (t) _ J 10 J Ql I 1 'o fo\ 

Po Jt)At) _At) At)' Pl At) At) _At) At)' 

'0 °10 "OWIO '1 u 01 a 10J01 

N.B. If r^^o - b$f$ = or r?h$ - 6$/oi = ° 5 there exists a sin S le special 

• i (*) CO n 

point Po — Pi = 0. 

The first special point — p± — corresponds to the extinction of the 
(two) populations. We note that this special point, when reached after some 
trajectory depending on the initial conditions, does not move in phase space: it 
is a fixed point. The second special point is much more interesting. In general, 
the coefficients in the model equations ^ can change their values in the 
course of time. Thus, the second special point can move in the phase space. 
As the phase point corresponding to the system state moves, one can observe 
an interesting hunt for equilibrium: the phase point 'hunts' the special point. 
The dynamics of the system trajectory in the phase space can become quite 
complicated, as seen in Fig. 1, for a particular numerical example. Changing 
ri leads to additional movement of the special point as shown in Fig. 2 . 

As we can see the system trajectory (marked with the dashed line in Fig. 
1) has to make a (long) journey through the phase space before finding and 
reaching the non-zero state described by the moving special point. This occurs 
despite the fact that at the initial time the system state and the special point 
are not very far one from another. Thus, the possibility that a time change in 
system coefficients enriches the dynamics in comparison to the case of constant 
coefficients of the model equations is easily demonstrated. 

Quite interesting are the phenomena connected to the linear stability of the 
two special points from Eqs.([7|-([8]). The eigenvalues connected to the stability 
of the special point (0, 0), which is a fixed point, are 



At) _ 1 
Ai, 2 - g 



^+rf ,) ±V(rS°-r? ) ) a + 4/«/S> 



(9) 



Let us assume that the coefficients connected to the mass media conversion are 
both positive, i.e. /q*' > and f[g > 0. Then the two eigenvalues are real. If 
the system parameters are constant in time and | Ai |< 1 and | A2 |< 1 the 
special point is stable. But as the parameters can vary with time the condition 
I Ai |< 1 and | A2 |< 1 can be violated and the point may become unstable for 
some time interval. If it is unstable, it cannot be reached. 

In addition a special point may move periodically; whence any attempt for 
the system trajectory to catch the moving special point can lead to a complex 
motion of the phase trajectory in the phase plane. This can be illustrated on 
the basis of the following example. Let / 01 = and 610 = and fix the values 
of parameters as mentioned in the caption of Figs. 3-4. Let all parameters 
be constant, except for ro which oscillates. The result is an oscillating special 
point; the system trajectory moves in a complex way in the phase space (see 
Fig. 3). 
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Of course, there can be more than one ideology, and the so called non- 
believers, present in the region. Then the existing ideologies can compete for 
adepts, such that the dynamics of the state point in the phase space can become 
more complicated. An illustration of such a case is shown in Fig. 5 for a 
particular numerical example. 

3 The case of two regions 

Consider two regions, I and II, and some possibility for migration of believ- 
ers or non-believers between the regions. Let us note the n + 1 densities of 
populations/opinions in the first region as pi,o,Pi,i, ■■■,pi, n and the densities of 
populations in the second region as pn,o, Pn,i, • ■ • , Pu,n (some of p\^ and pnj 
can be 0). Let it be conventionally accepted from now on that the label before 
the comma in the index refers to the region, while the index after the comma 
refers to the ideology type or population. The leaving from and coming to a 
region are accounted for by introducing new terms in the model equations with 
appropriately defined specific diffusion coefficients. Let the coefficient of leaving 
(L) the area I toward the area II by adepts of the i-th. population be ii-j.ii,. and 
the coefficient of leaving of the area II toward the area I by the adepts of the 
i-th population be Ln-s-i, .- The corresponding coefficients for the incoming (/) 
followers of the i-th population are /n_>.i and /i_>h h , or ii<_ii,j and In^i,iFl A 
necessary balance exists between the number of adepts leaving one area when 
going into the other, implying some obvious relationships between these coef- 
ficients. We assume that an adept does not change his/her opinion by simply 
leaving from or going into a region. 

The model equations for the two area (n+1) population number evolutions 
become 

n n 

Jt+i) _ (t) _ (t) (t) At) (t) At) ft) <t) _ 

Pl,i P\,i — r l,iPl,i + h,ijPl,j + / j a l,ijPl,iPl,j 

3=0 j=0 

At) n (t), T (t) (t) nm 

L I-+lI,iPl,i + h^-ll.iPll.i \ w ) 

n n 

„(*+!)_ (t) (t) (t) ,y- f M (t) ,y t W At) (t) _ 

3=0 3=0 

r(t) Jt) , At) At) ni s 

Let us consider, for illustrating the notations, the particular case of exchange 
between two areas when only two populations/opinions are present in each re- 
gion: the non-believers po and followers of the (so called prevailing) ideology 1, 

3 It might be interesting to have one or several indices attached to L and / in order to 
distinguish among the migration causes like migration currents due to external forces or eco- 
nomic reasons due to gradients in standard of living, though this would make the notations 
rather elaborate. Such a remark is highly appropriate for the cases discussed below. We thank 
one of the reviewers for such a good suggestion. 
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i.e. pi. The model equations are reduced to 



_ Jt) _ Jt) Jt) , ,(*) (t) , ,(t) (t) (t) _ At) it) At) (t) 
n,o fi,o ~ 'i.oPifl < J 1,01^1,1 ' "1,01^1,0^1,1 ^1-^11,0^1,0 ' - t i«-n,orii,o 

At+i) _ Jt) _ Jt) Jt) M) Jt) Jt) Jt) Jt) _ At) Jt) Jt) Jt) 
fn.,1 fi,i ~ 'i, 1^1,1 t J i, 10^1,0 "r "1,10^1,1^1,0 ^1^11,1^1,1 ' J i<-ii,irn,i 

-(*+!)_.{*) _J*)J*) + /(*) J.*) 4.i,(*) J*) J*) _r(*) „(*) 4.rW „(*) 

^11,0 /''II.O — 'II,0PII,0 ' J II, 01^11,1 ^ "II, 01^11,0^11,1 ^11^-1,0^11,0 ' J 11^1,0^1,0 

„(*+!)_„(*) _ J*) J*) .A*) Jt) , h (t) Jt) Jt) _r(t) Jt) At) Jt) 
PII,1 VlIA — '11,1^11,1 ^ J II, 10^11,0 ^ "II, 10^11,1^11,0 -^11^1,1^11,1 "r i II^I,lFI,l 



(12) 



From the numerous possible scenarii resulting from the above, we shall fur- 
ther restrict the present considerations to two cases only: 



massive leaving of one region by non-believers (i — 0) under the pressure 
of the followers of the prevailing ideology (i = 1); e.g., protestants leaving 
Europe because of catholics pressure in the XVI -th century or so. 

massive invasion of one region by the followers of the prevailing ideology 
(i = 1) from the other region; e.g., catholics conquistadors going to South- 
America in the XV-th century or so. 

N.B. Although both cases appear at first rather equivalent, the two cases under 
interest should nevertheless be suggestively markedly different to the astute 
reader, because of the birth/death rate conditions imposed on populations 
and 1 respectively , from the mere beginning. 



3.1 Mass departure 



Let the followers of the prevailing ideology (i = 1) put pressure on the non- 
believers (i = 0) in order to fully convert some region toward their ideology; 
thus let the a massive leaving of one region, e.g. I, by non-believers (i = 0) 
toward the other region, i.e. II. We assume that no exchange of populations 
between the two regions happens except for the leaving of non-believers from 
region I to region II, i.e. and all other coefficients L and /, except 



of course Jj q, in Eq. (|12), equal to 



From the first equation 



the number of non-believers (i — 0), in I, is 



of the system (12 1, the condition for a decrease in 



At) 

-^i->ii,o 



> r [t) + 

> 'i,o + 



(*) 

At) 

h,oi ( t ) 



b {t) o (t) 

°1, 01^1,1' 



(13) 



Eq. ( 13 1 tells us that if the coefficient describing the leaving is large 

enougnin order to compensate the birth rate and the effect of the conversion 
mechanisms, then the population of non-believers in the region I will start to 
decrease. After some time, the followers of the ideology (i = 1) can thus become 
the majority of the population in that region. 
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Due to the balance condition, i.e. the conservation of the number of non- 
believers in both regions, the number of non-believers in the region II will in- 



crease. Thus, from the third equation in (12), one has similar condition on 

r(t) 

^i-m.i 
II, i.e. 

,W Jt) _L fit) At) j_ At) jt) Jt) 

hi.: 

(t) 



Lj\jj , but relevant for the behavior of populations (i — and i = 1) in region 



T (t) T (t) '11,0^11.0 WII.Ol^II.l ' 01^11.0^11,1 ( -, A x 

^1-5-11,0 = J II^I,0 > (t) • ! 14 ' 



3.2 Mass invasion 

Another scenario leading to some domination of a region by an ideology may 
be due to the invasion of that region from another region by the followers of 
the prevailing ideology in the latter region. Notice that this scenario is more 
complex than the scenario of leaving in Sec. |3.1| as it requires common efforts 



by the followers of the ideology in the two considered regions. Let, in Eq. ( 12 ), 
Ij^u 1 = i 7^ and all other coefficients L and I in jl2| be 0, except of 

course Ljf^ v The condition for an increase in the number of believers (i = 1), 
in region I, is 

Jt) At) At) (t) At) (t) (t) 

T (t) 1,1 "l,l ' JI,10"l,0 ' l,10"l,l"l,0 / 1t -x 

J i<— n,i > <t) • y i0 ) 

Pu'l 

In other words, the invasion has to compensate for the death rate (remember 
that r\ -{) is negative) and is supported by the conversion mechanisms, through 
/j j and bfl Q . Thus, if the followers of an ideology manage to organise mass 
invasion combined with massive media propaganda (large fj^) and large con- 
versation conversion (large 6j 10 ) an d if they face a "tired" population of non- 
believers (small or even negative r[*g ) which does not take measures to decrease 

the number of followers of the ideology by conversion (small f^L and b\ q-J then 
the invasion could be successful. The balance condition, i.e. the conservation of 
the number of ideology (i = 1) adepts in both regions, implies that the number 



of the latter in region I will increase. Thus, from the fourth equation in ( 12 1 , one 
has similar condition on L 1 ^^ 1 , 
(i = and i = 1) in region II, i.e. 



has similar condition on 1 , but relevant for the behavior of populations 



(*) 

rW > Jt) , At) Pll,0 ■ At) it) , lfi v 

-^11— 5-1,1 ' 11,1 ' ./II, 10 (t) ^ "11,10^11.0 • 

Piii 
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4 Arbitrary number of regions. Further com- 
ments 



The final complication of the model with respect to the number of discussed 
regions is to consider arbitrary number of regions. We separate the entire studied 
two-dimensional area into Sk regions, k = I, II, . . . , K. In each region we assume 
i = 0, 1, ... ,n populations with densities pk.i where i = denotes the population 
of non-believers in any ideology. In general we allow migration of members of 
populations between any two of the areas even between those that do not have 
common boundary. The corresponding exchange coefficients are denoted as 
Lk^fi.i an d Ik<r-l,i- The system of K ■ {n + 1) model equations read^] 

n n 

Jf+i) _ J.t) _ (*)-(*) + V f (t) o (t) + V b {t) o (t) o (t) ~ 

3=0 3=0 

K K 

^ ^k-yl.iPu + ^ h<-LiPlA 
1=1 1=1 

i = 0,l,...,n; k = 1,11,..., K (17) 

This system describes the dynamics in a K • (n + 1) - dimensional phase space 
where numerous special points and complicated attracting manifolds arise, van- 
ish, and move with time. The system trajectory can be very complicated as 
it travels in a high-dimensional phase space in order to come close to some 
attracting manifold. 



5 Practical Remarks 

At this stage, many interesting research questions arise about applications: for 
example, questions on the control of the dynamics of the ideological interaction. 
Let us discuss a situation on the control of the spreading of ideologies in a 
selected region of a territory consisting of several sub-regions on the basis of the 



general model system Eqs.(17). For bettering such a control, the dimension of 
the phase space of the system should be reduced. This can happen by regulating 
and reducing the inflow of followers of unwanted ideologies, i.e. technically, by 
lowering the values of I^Lu { - An outflow of followers of unwanted ideologies can 
be stimulated, if this does not affect the coefficients of mass media conversions. 
Otherwise the outflow has to be reduced also, i.e. lowering Lu\i v Without 
significant inflow and outflow in the selected region, the dimension of the phase 
space of the ideological dynamics is practically reduced. Within a purpose of 
control, one has to deal with the followers of the ideologies and with the non- 
believers in the region. 



4 It is surely obvious to the reader that a few coefficients are necessarily equal to 0, but all 
terms are hereby written in order to maintain the number of notations at a reasonable level. 
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The dynamics of the corresponding populations is determined by the birth 
and death rate of the non-believers as well as by the death rates of the followers of 
all the other ideologies and by the conversion coefficients. Anyone's experience 
shows that the conversion by conversation is much part of an e.g. religion 
conversion in a dynamical context [23 (Ml HS1 130] , while the conversion 

due to media impact is well known to be the primary cause of evolution in 
politics [311 E2 E3 EH] ■ Of interest is also the conversion of an opinion about 
goods and consumer products [3_3 ESI EH [38]. For completeness, recall that 
external (propaganda) influence and socio-economic pressure can be mapped 
into social forces |39j . The correspondence with our approach is left for another 
publication. Nevertheless, they would be easily apparent if one was developing 
the rhs of the complete dynamical equations 

(t+2) _ (t+1) (t) 
Pk,i ^Pk,i + Pka — • • • i 

i = 0,l,...,n; fe = I,H, ...,K (18) 

Thereafter, one can discuss qualitatively cases of recent ideological control, for 
example in Arab countries and in China. Massive invasion will be considered 
at an earlier time, i.e. the XlV-th century invasion of Europe by Turks. 

Recently, the failure of ideological control has been seen in several Arab 
countries. Several reasons for this come in mind on the basis of the discussed 
model. Let in this case the index i in the model equations denotes: i — - no 
ideology; i — 1 - official ideology; i = 2 - Islamism; i — 3 - liberal ideology; 
i = 4 - moderate Islam. The official ideology connected to corresponding dicta- 
torship regimes has gained followers mostly through the domestic mass media 
propaganda and conversations supported by financial support)^] This leaded to 
significant positive value of /io- &io had significant positive value too. Islamists 
were not successfully converted to the official ideology, i.e., the values of /i 2 
and &12 were small and negative. More successful was the conversion from the 
liberal to official ideology which means that /13 and 613 had positive values. 
Relatively successful was the conversion from Islam to the official ideology (as 
the official ideology was partially based on Islam). As consequence of this the 
coefficients /14 and 614 had significant positive values. Thus the strategy of the 
corresponding government was to convert as many as possible non-believers and 
followers of the Islam to the official ideology and to decrease as much as possible 
the conversion of the followers of the official ideology to Islamism. As well as 
the government had enough money and propaganda power it was possible to 
follow such a strategy. 

The liberal ideology spread by the mass media was mainly connected to the 
west technology: Interned and West TV channels and military shows. This 

5 There is also the "service and welfare support" to the need of the population, - which is 
not really mass media propaganda nor conversation, but rather an opinion "formation" due 
to some "theoretically gratuitous altruism" 40 41, 42, 43 44 by "charitable organisations" 
following resources transfer, with political goals 45 46, 47, 48 49 . In view of its clustering 
effect, somewhat a group response to group action, it might be considered to be a higher order 
term in Eqs.lpl, or still be integrated into each f?p coefficient. 
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ideology was not able to compete successfully against the government ideology 
which means that the corresponding conversion coefficients /3o,&3o had smaller 
values in comparison to /i , 610 • The conversion from official ideology to liberal 
ideology was not successful, i.e., /31 and 631 had negative values. The values 
of /32 and 632 were very small negative (there was no conversion of Islamists 
to liberal ideology and small number of followers of liberal ideology converted 
to Islamism). The conversion of the followers of the Islam to liberal ideology 
was also not successful. Thus the values of / 3 4 and 634 were almost zero. Thus 
the strategy of the followers of the liberal ideology was to convert the non- 
believers (mainly young people) and to minimize the decreasing of the followers 
by conversion to the official ideology and Islamism. 

The Islamism gained followers by conversations, hypocritical altruism and 
the domestic mass media. Using the power of the Islam large values of /20, /24 
and 620^24 have been maintained. /2i,&2i had small positive values as some 
supporters of the official ideology converted to Islamism. The values of / 2 3 and 
623 were also positive but much smaller as very few supporters of the liberal 
ideology converted to Islamism. 

The followers of the moderate Islam increased in the course of years because 
of the large conversion rates /4o,&4o from (youth) non-believers to followers of 
the Islam. Much smaller conversion happened from Islam to followers of the 
official ideology (small negative values of /41 and 641) and even smaller was 
the conversion from the followers of the liberal ideology (very small positive 
values of and 643). The conversion of fanatic Islamists to the moderate 
Islam was not successful (negative values of /42 and 642). The centuries old 
strategy of the followers of moderate Islam to convert very successfully (young) 
non-believers worked very well and compensated the losses due to conversion to 
official ideology and Islamism. 

The outflow of people from these countries toward others was not (at this 
time) large enough, though the large birth rates in the Arab countries produced 
a large amount of jobless youth which becomes resistive to accept the official 
ideology in the corresponding Arab country. Conversion to the corresponding 
official ideology by conversation, was apparently not very successful either and 
the values of /10 and 610 remained positive but decreased very much with re- 
spect to the values of the corresponding coefficients for the other ideologies. 
Thus in the course of the years the non-believers were converted much more 
successfully to moderate Islam, Islamism and even to the liberal ideology. The 
spreading of liberal opinions following mass media conversion, such as Internet, 
Facebook, Twitter, etc., was not totally controlled and thus the non-believers 
were even more resistive for conversion to the official ideology. Thus the official 
(and mostly totalitarian) ideology in these countries had decreased percentage 
of followers and finally collapsed. Now in these countries we observe fast growth 
of the fanatic Islamists, many moderate followers of Islam, still influential popu- 
lation of followers of the former official ideology and a slowly increasing amount 
of people that "convert" to some liberal ideology, under the influence of the 
mass media from the West. Then the future development in these societies in 
the next years will depend much on the competition between the liberal ideol- 
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ogy, Islamists and moderate Islam followers. The liberal ideology has a problem: 
it can not convert very successfully the followers of the moderate Islam. The 
Islamists do not have this problem. Thus a possible capsulation of the Arab 
societies (and especially the decreasing influence of the West mass media) could 
lead to a dominance of Islamists and Islamic ideology which is quite deeply 
rooted in these societies and has a large tradition of conversion of non-believers. 

Much more successful is the ideological control in China. Many elements of 
the capitalist economic system help for successfully absorbing the inflow of young 
people who have studied in the West. This absorption is easily transformed 
by (social or economic or political) pressure conversion in the official Chinese 
ideological system. Birth rates are much under control and the outflows of 
followers of different ideologies are much regulated [5CT . 

Let in the case of China the value of the index i denote ideologies as follows: 
i = - no ideology; i — 1 - official ideology; i = 2 - liberal ideology. Mass media 
is under control and the Chinese government takes steps toward controlling 
the Internet and the connections between social networks. This ensures large 
positive values of the conversion coefficients /io and &io whereas the values of 
the corresponding coefficients for the liberal ideology /20 and &20 are positive 
but much smaller. As the official ideology contains element of liberalism the 
values of conversion coefficients /12 and 612 are small positive and the values of 
corresponding coefficients for the liberal ideology / 2 i and 6 2 i are negative. We 
note that in principle the conversation conversion is much out of controj^] but 
up to now it is successfully countered by governmental propaganda and internal 
social pressure and the value of 620 is small. Thus the ideological situation 
in China is stable and the ideological dynamics is dominated by the official 
ideology. As no large demographic problem is expected in the next decades and 
if the economic development continues, the stability of the ideological situation 
is China is expected to remain largely unchanged. 

Finally, from the historic situations connection to mass invasion of follow- 
ers of some ideology to a region dominated by another ideology let us consider 
the invasion of the Ottoman Turks into Europe that started from the second 
half of the 14-th century. We consider two regions: region I.: the Balkans (the 
invaded region); and region II. (Anatolia, Asia Minor and surrounding areas 
of Asia) . At this time the invaded region: the Balkans were populated by fol- 
lowers of Christianity (mostly from the orthodox branch but with presence of 
followers of Catholicism too). Weakened by the 'Black death' plague the small 
Balkan countries were conquered one after another by the Muslim Ottoman 
turks who started to impose Islam, ruining churches or transforming these into 
mosques. Let us discuss qualitatively the situation from the point of view of 
the model presented above. Let the values of the index i mean the following: 
i = - non-believers; i = 1 - Christians; i = 2 - followers of the Islam. The 
birth rates supported the population of non-believers that in course of years 
converted to Christianity or Islam. At that time there was no mass media as 

6 Some investigation in such a respect is worth to be pointed out: in the Moss model 51 , 
social embededness takes the form of observation of neighbours public consumption activities 
such as garden watering and car washing as well as word of mouth communication. 
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we know it today so /i,oi, /n,oi, /i,icb /n,io were practically and there was 
no possibility of media influence on non-believers in order to convert them to 
Christianity. The values of the coefficients /1.02, /n,02, /1.12, /ii,12 were practi- 
cally too as there was no conversion of the followers of Islam to Christianity or 
to the population of non-believers. The public spreading of the administrative 
orders of the Ottoman administration and the lower taxes for Muslims can be 
treated as sort of media that contributed to conversion and led to small positive 
values of the coefficients /i,20i /n,20>/i,2i and /n,2i- The main process for en- 
suring the existence of the Christian population was the conversion of (young) 
non-believers to Christianity by means of conversation, i.e., the value of 6t,io 
was large positive. The value of 611,10 was small positive as in the region II 
the Christian population was surrounded by large Muslim population that put 
pressure acting in the direction of muslimisation of the non-believers in that re- 
gion. This was a branch of the Ottoman assimilation mechanism that included 
displacement from Europe to Asia and then muslimisation. In the two regions 
there was almost no conversation conversion from Christianity or Islam to the 
population of non-believers. Because of this the values of the coefficients 61,01 > 
bi,02, bu,oi, 611,02 were negligibly small. The conversation conversion from Islam 
to Christianity was also negligibly small, i.e., the values of the coefficients 61,12, 
&ii,12 were practically 0. 

The politically dominating Muslim population was not able to keep large 
enough level of conversation conversion in the region I: the values of 61,20 and 
61,21 were positive but far away from the needed values for fast conversion of 
the Christian population to Islam. The population of entire regions refused to 
convert to Islam. For increasing the value of /i,20 the 'enichars'' tactic was 
used for almost three centuries - Christian children were taken away from their 
parents and then converted to Islam. The increase of the values of /i,2i and 
61,21 was attempted by campaigns to muslimise entire regions of the Balkan 
peninsula (todays Albania where much more than the half of population were 
catholics, todays Bosnia and Herzegovina, and todays Rodopa mountain region 
in Bulgaria where the population consisted of almost 100 % followers of Ortho- 
dox Christianity). Who refused to convert was killed. In the region II (Asia) 
the situation was different as there the Muslim population was large majority 
and the pressure for conversion of Christian population was effective, i.e., the 
values of 611,20 and 611,21 were positive and large enough for the conversion to 
happen. 

Important point of the Ottoman strategy for muslimisation of the Balkans 
was connected to migration. There was massive migration of Muslims from 
region II (Asia) to region I (Balkans) especially after Timur (known also as 
Tamerlane) defeated and captured the Ottoman sultan Bayezit I in 1402 near 
todays Ankara. The remains of the Ottoman army were transported from Ana- 
tolia (region II) to Thrace (region I) by Venetian shipsQ Further migration 

7 Venetians preferred to have the weak (at that time) Ottoman turks as enemy and a buffer 
against the mighty Timur but Timur aimed to restore the Mongolian Empire of Chengis Khan 
and turned to conquer China. The Ottoman army could reorganize and crushed the Bulgarian 
uprising in 1405. Thus the way to Europe was opened for the Ottoman that managed to 
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followed and the coefficients Ln-»i,2 and 11,2 had large positive values. The 
migration of Muslims in opposite direction was very small and the values of the 
coefficients L1-+112 and im-1,2 were practically 0. 

Parallel to the migration of the Muslim population the Ottoman turks dis- 
placed the Christian population and non-believers (children) and part of it was 
sent from region I to region II (example of this is the displacing of large part of 
Christian population of Belgrade to Istanbul) so the coefficients Li-»ri.i, Li-sj^o, 
hi^i.a , hi^i.i had also large positive value. Such population usually remained 
in Anatolia and had been converted to Islam through the centuries of Ottoman 
rule. Thus the values of the coefficients Lt^t 1, Ln^i.o, hi-iifi and h^n,i were 
very small. 

The conversion strategy of the Ottoman turks acted for a very long time. 
Because of this, 450 years after the start of the invasion between 30% and 50% of 
the population of several regions of the Balkans has been practically converted 
to Islam. 



6 Conclusions 

In this paper we have discussed a discrete model of ideological interaction which 
is easy for practical implementation for prognosis of the dynamical evolution of 
ideological situations. The model system is based on a discrete time process 
within a Lotka-Volterra framework, but allowing for migration of populations 
between areas. The model contains specific terms mimicking as much as possi- 
ble stylised facts or, as usually, common sense in opinion evolution dynamics. 
Each country or region of countries can be separated into a number of regions 
for which some data on the number of followers of different ideologies can be col- 
lected. Birth and death rates are often available from statistical sources indeed. 
From these sources estimations of the coefficients of mass leaving and incoming 
can be usually obtained as well. What remains, and is the most difficult task, 
is the determination of the conversion coefficients. However it seems that they 
can be obtained by standard or appropriate statistical procedures, like opinion 
pool evaluations, or by indirect means of behaviour observation. Of course long 
time interval observations are required in order to determine the time change of 
the coefficients for model equations. 

After determining the coefficients one can launch the model, starting from 



Eqs.(17), with known or appropriate available initial conditions. Under the 
assumption that the coefficients of the model equations vary rather slowly, it 
can be fairly assumed that either they are constant during several time steps or 
they change following a linear law during these steps. Therefore, the model is 
expected to give a valuable approximate picture of the ideological situation in 
the observed region for these time steps. Of course, such approximations can be 
best applied when the situation in the corresponding regions is stable and when 
no large social change happens. In times of large social changes, and extreme 



besiege even Vienna two times in the next 300 years. 
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events, like wars, revolutions, the constant or linear approximation can be 
unsatisfactory. 

As discussed here above, the model admits a very large number of possible 
scenarios. In conclusion, even without solving in detail the system of equations, 
one can obtain some orientation about ideology situations in different regions 
merely on the basis of qualitative discussions of the equations of the model. 
Only limited to opinions, it can be fruitfully compared with attempts of the so 
called World Model [52] for socio-economic considerations. 

Acknowledgments 

We would like to acknowledge the stimulating discussions with many colleagues 
participating in the COST MP0801 Action 'Physics of competition and conflict'. 
We thank the COST MP0801 Action for the support of our research. We thank 
the reviewers for questions allowing to improve the manuscript. 



References 

[1] Vitanov, N.K, Dimitrova Z.I., and Ausloos M., Verhulst-Lotka-Volterra 
(VLV) model of ideological struggle. Physica A 389 (2010) 4970 - 4980. 

[2] Hinich, M.J. and Munger M.C., Ideology and the Theory of Political Choice 
(The University of Michigan Press, Ann Arbor, 1996). 

[3] Montroll, E.W. and Badger W.W., Introduction to Quantitative Aspects of 
Social Phenomena (Gordon and Breach, New York, 1974). 

[4] Brauer, F. and Castillo-Chavez, C, Mathematical Models in Population 
Biology (Springer, Berlin, 2001). 

[5] Vitanov, N.K., Jordanov, LP., and Dimitrova, Z.I., On nonlinear dynamics 
of interacting populations: Coupled kink waves in a system of two popu- 
lations. Communications in Nonlinear Science and Numerical Simulation 
14 (2009) 2379 - 2388. 

[6] Hastings, A., Population Biology: Concepts and Models (Springer, New 
York, 1997). 

[7] Vitanov, N.K., Jordanov, LP., and Dimitrova, Z.I., On nonlinear popu- 
lation waves. Applied Mathematics and Computation 215 (2009) 2950 - 
2964. 

[8] Okubo, A., and Levin, S., Diffusion and Ecological problems: Modern Per- 
spectives (Springer, New York, 2001) 

[9] Vitanov, N.K., Dimitrova, Z.I., and Kantz, LL, On the trap of extinction 
and its elimination. Physics Letters A 349 (2006) 350 - 355. 



16 



[10] Kaplan, D. and Glass, L., Understanding Nonlinear Dynamics (Springer, 
New York, 1995). 

[11] Vitanov, N.K., Breather and soliton wave families for the sine-Gordon equa- 
tion. Proceedings of the Royal Society of London A 454 (1998) 2409 - 2423. 

[12] Mees, A. I., Nonlinear Dynamics and Statistics (Birhkauscr, Boston, 2001). 

[13] Vitanov, N.K. and Yankulova, E., Multifractal analysis of the long-range 
correlations in the cardiac dynamics of Drosophila melanogaster. Chaos, 
Solitons & Fractals 28 (2006) 768-775. 

[14] Panchcv, S., Spassova, T., and Vitanov, N.K., Analytical and numerical 
investigation of two families of Lorenz-like dynamical systems. Chaos, Soli- 
tons & Fractals 33 (2007) 1658-1671. 

[15] Kiel, L.D. and Elliott, E. (Eds.)., Chaos Theory in the Social Sciences: 
Foundations and Applications (The University of Michigan Press, Ann Ar- 
bor, 1996). 

[16] Dimitrova, Z.I. and Vitanov, N.K., Dynamical consequences of adapta- 
tion of the growth rates in a system of competing populations. Journal of 
Physics A: Mathematical and General 34 (2001) 7459-7473. 

[17] Boeck, T. and Vitanov, N.K., Low-dimensional chaos in zero-Prandtl num- 
ber Benard-Marangoni convection. Physical Review E 63 (2002) 037203. 

[18] Dimitrova, Z.I. and Vitanov, N.K., Chaotic pairwise competition. Theoret- 
ical Population Biology 66 (2004) 1-12. 

[19] Gardiner, G.W., Handbook of Stochastic Methods for Physics, Chemistry, 
and the Natural Sciences (Springer, Berlin, 1983). 

[20] van Kampcn, N.G., Stochastic Processes in Physics and Chemistry (Else- 
vier, Amsterdam, 2007). 

[21] Kanz, H., Holstein, D., Ragwitz, M., and Vitanov, N K., Markov chain 
model for turbulent wind speed data. Physica A 342 (2004) 315-321. 

[22] Vitanov, N.K., Upper bounds on the heat transport in a porous layer. 
Physica D 136 (2000) 322-339. 

[23] Vitanov, N.K., Convective heat transport in a fluid layer of infinite Prandtl 
number: upper bounds for the case of rigid lower boundary and stress-free 
upper boundary. European Physical Journal B 15 (2000) 349-355. 

[24] Vitanov, N.K., Upper bounds on convective heat transport in a rotating 
fluid layer of infinite Prandtl number: case of intermediate Taylor numbers. 
Physical Review E 62, (2000) 3581-3591. 



17 



[25] Hayward, J., Mathematical modeling of church growth, Journal of Mathe- 
matical Sociology 23 (1999) 255-292. 

[26] Hayward, J., A general model of church growth and decline, Journal of 
Mathematical Sociology 29 (2005) 177-207. 

[27] Ausloos, M. and Petroni, F., Statistical dynamics of religions and adher- 
ents, Europhysics Letters 77 (2007) 38002. 

[28] Ausloos, M. and Petroni, F., Statistical dynamics of religions, Physica A 
388 (2009) 4438-4444. 

[29] Ausloos, M. and Petroni, F., On World religion adherence distribution 
evolution, in Takayasu M., Watanabe T., Takayasu H. (eds.) Econophysics 
Approaches to Large-Scale Business Data and Financial Crisis (Springer, 
Berlin, 2010) pp. 289-312. 

[30] Ausloos, M., On religion and language evolutions seen through mathemat- 
ical and agent based models, in Rangacharyulu, C. and Haven, E., (Eds.) 
Proceedings of the First Interdisciplinary CHESS Interactions Conference 
(World Scientific, Singapore, 2010) pp. 157-182 

[31] Converse, P.E., Changing conceptions of public opinion in the political 
process, The Public Opinion Quarterly 51 (1987) S12-S24. 

[32] Iyengar, S., Is Anyone Responsible?: How Television Frames Political Is- 
sues (University of Chicago Press, Chicago, 1994) 

[33] McQuail, D., The influence and effects of mass media, in Curran, J., Gure- 
vitch M., and Woolacott J. (Eds,) Mass Communication and Society (Sage 
Publications, London, 1977) pp. 70-93. 

[34] Gerber, A.S., Karlan, D., and Bcrgan, D., Does the media matter? A 
field experiment measuring the effect of newspapers on voting behavior 
and political opinions, American Economic Journal: Applied Economics 1 
(2009) 35-52. 

[35] Turnbull, P.W. and Meenaghan, A., Diffusion of innovation and opinion 
leadership, European Journal of Marketing 14 (1980) 3-33. 

[36] Cheung, C.M.K., Lee, M.K.O., and Rabjohn, N., The impact of electronic 
word-of-mouth: The adoption of online opinions in online customer com- 
munities, Internet Research 18 (2008) 229-247. 

[37] Buttle, F.A., Word of mouth : understanding and managing referral mar- 
keting, Journal of Strategic Marketing 6 (1998) 241-254. 

[38] Bass, F.M., A new product growth model for consumer durables, innovative 
and imitative behavior. Management Science 18 (1969) 215-227. 



18 



[39] Montroll, E.W., Social dynamics and the quantifying of social forces, Pro- 
ceedings of the National Academy of Sciences USA 75 (1978) 4633-4637. 

[40] Frolich, N., Oppenhcimer, J., and Young, O., Political Leadership and Col- 
lective Goods (Princeton Univ. Press, Princeton, 1971). 

[41] Krebs, D.L., Altruism- an examination of the concept and review of the 
literature, Psychological Bulletin 73 (1970) 258-302. 

[42] Underwood, B. and Moore, B., Perspective-taking and altruism Psycholog- 
ical Bulletin 91 (1982) 143-173. 

[43] Gintis, H., Bowles, S., Boyd, R., and Fehr, E., Explaining altruistic behav- 
ior in humans, Evolution and Human Behavior 24 (2003) 153-172. 

[44] Fehr E. and Fishbacher U., The nature of human altruism, Nature 425 
(2003) 785-791. 

[45] Kuran, T., On the notion of economic justice in contemporary islamic 
thought International Journal of Middle East Studies 21 (1989) 171-191. 

[46] Ayubi, Nazih N.M., Political Islam: religion and politics in the Arab world 
(Routlcdge, London, 1991) 

[47] Giilalp, H., Globalization and political Islam: The social bases of Turkey's 
Welfare Party, International Journal of Middle East Studies 33 (2001) 433- 
448. 

[48] Bayat, A., Activism and social development in the Middle East Interna- 
tional Journal of Middle East Studies 34 (2002) 1-28. 

[49] Roscfsky Wickham, C, Mobilizing Islam: Religion, Activism and Political 
Change in Egypt (Columbia Univ. Press, New York, 2002) 

[50] Yuan Tien, H., China: Demographic Billionaire. Population Bulletin 38 
(1983) 1-48. 

[51] Downing, T.E, Moss, S. and Pahl Wostl, C, Understanding climate policy 
using participatory agent-based social simulation, in Moss S. and Davidsson 
P. (Eds.) Multi Agent Based Social Simulation (Springer, Berlin, 2001) pp. 
127-140. 

[52] Meadows, D.H., Meadows, D.L., Randers, J., and Behrens III, W.W., The 
Limits to Growth, A Report for THE CLUB OF ROME'S Project on the 
Predicament of Mankind (Universe Books, New York, 1972). 



19 



400 
300 
P^OO 
100 


-50 50 100 150 200 

Po 



Figure 1: Illustration of the system trajectory and the special points in the 
phase space for the case of a single area and two populations, - non-believers and 
followers of an ideology. Solid line: movement of the non-zero special point given 
by Eq.(|8|. Dashed line: system trajectory calculated on the basis of Eq.([5]) and 
Eq.([6]). The parameters are appropriately chosen in order to illustrate how the 
system phase point 'hunts' and 'catches' the (non trivial) special point {pq\ pi) 
as t — > oo. Initial values: po = 0.98, p\ = 0.02. The constant parameters are: 
/$ = 0.0002; = 0.0001; rf 5 = -0.001; 6$ = -0.0001; 6$ = 0.0001. All 
such parameters are assumed to be constant. The only time changing parameter 
is: = 0.01 tanh(O.Oli). The maximum value of t used is t = 40000. 
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Figure 2: Effect of the change of r^. Solid line: the corresponding trajec- 
tory from Fig.l. Dashed line: the system trajectory for = — 0.01[1 — 
0.01 tanh(O.OOli)]. Illustration of the system trajectory and the special points in 
the phase space for the case of a single area and two populations, - non-believers 
and followers of an ideology. Solid line: movement of the non-zero special point 
given by Eq.(|8]). Dashed line: system trajectory calculated on the basis of 
Eq.([5]) and Eq.([6]). The parameters are appropriately chosen in order to illus- 
trate how the system phase point 'hunts' and 'catches' the (non trivial) special 
point (fig ,pi ) as t — > oo. Initial values: po — 0.98, pi — 0.02. The constant 
parameters are: f$ = 0.0002; = 0.0001; 6$ = -0.0001; 6$ = 0.0001. All 
such parameters are assumed to be constant. The only time changing parameter 
is: = O.Oltanh(O.Olt). The maximum value of t is * = 40000. 
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Figure 3: Illustration of a periodic motion of the non-trivial special point, 
i.e., motion of the components p$\ p^ of the special point from Eq. |sj. Solid 
line: p^; dashed line: pf> . Initial values: p^ = 6.0 • 1CT 2 ; p { ? ] = 2.0 ■ 10" 1 ; 
the parameters /io = 10~ 4 ; /oi = 0; r\ = — 10~ 3 ; &oi = —0.1; feio = are 
assumed to be constant, but r = 1.5 tanh(10 _3 i) — 1.05 sin(10 _3 i). N.B. the 
y-axis should read po,i and not p\p. 
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Figure 4: Illustration of the corresponding complicated trajectory of the system 
in the phase space from Eq. All values of the parameters and quantities 

are chosen in order to make the illustration more visible. Initial values: p °^ = 
6.0 • 10~ 2 ; pf> = 2.0 ■ 10 -1 ; the parameters f w = 10~ 4 ; f 01 = 0; n = -10~ 3 ; 
boi = —0.1; bio = are assumed to be constant, but r = 1.5 tanh(10~ 3 £) — 
1.05sin(10- 3 t) 
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Figure 5: Illustration of the trajectory of a system made of three populations, i.e. 
non-believers and followers of two distinct ideologies. The initial values for the 
run are: po = 1; pi = 0.2; p 2 = 0.1. Other parameters have the constant values: 
ho = 10- 2 ;/i2 = 10- 5 ;./oi = 0;/ 02 = 10- 6 ;/ 20 = 10- 3 ;/ 21 = lO" 3 ; ; n = 
-10" 3 ; r 2 = -lO- 3 ; b i = -0.1; b w = 0; b Q2 = -I0~ 4 ;b 12 - 10- 6 ;6 20 = 
10- 5 ;fe 2 i = -10~ 4 , but r = 1.5 tanh(O.OOlt) - 1.05 sin(0.001t) + 0.02 sin(0.03t). 
The trajectory is calculated up to t = 40000 
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